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In this article we provide a solution to Problem 2 in Gneiting, 2013. Specif- 
ically, we show that the even- resp. odd-dimensional Schoenberg coefficients 
in Gegenbauer expansions of isotropic positive definite functions on the sphere 
S d can be expressed as linear combinations of Fourier resp. Legendre coeffi- 
cients, and we give closed form expressions for the coefficients involved in these 
expansions. 

1 Introduction 



For an integer d > 1 we denote the unit sphere in Euclidean space R d+1 equipped with the 
! Euclidean norm by § rf := {x € R rf+1 : \\x\\ = 1}. A function h : S d x E> d — >■ R is positive 

' definite if 

in ■ 
oo 



2^ 2^ aidjhixiiXj) > 0, (1) 
i=\ j=i 



for all integers n > 1 and for every choice of constants a\ , . . . , a n el and every choice of 
pairwise distinct points x±, . . . , x n £ § d . The function h is said to be isotropic if there exists 
a function ib : \0, ir] — > R such that 

h(x,y)=if>(9(x,y)) 1 x,yeS d , 

where 6(x,y) = arccos((x, y)) denotes the great circle distance between x and y; here (.,.} 
is the standard scalar product in R rf+1 . 

We refer to *&d,d = 1,2,..., as the class of continuous functions tp : [0, ir] — s> R with 
i/)(0) = 1 for which the associated isotropic function h(x,y) = ip(8(x,y)) is positive definite. 
Therefore, ^ corresponds to the class of correlation functions for mean-square continuous, 



'The author's research was supported by the German Research Foundation (DFG) within the programme 
"Spatio/Temporal Probabilistic Graphical Models and Applications in Image Analysis", grant GRK 1653. 



1 



stationary and isotropic random fields on S rf (Jones, 1963). In addition, \P~t is denned as 
the class of elements if) G ^ for which the associated function h(x,y) is strictly positive 
definite, that is the sum in (1) is strictly positive. 

Isotropic positive definite functions on spheres have attracted research interest in several 
areas. They occur as correlation functions for stationary and isotropic random fields on the 
sphere and, hence, are studied in the context of spatial statistics (Banerjee, 2005, Huang, 
Zhang, and Robeson, 2011 or Hansen, Thorarinsdottir, and Gneiting, 2011). Furthermore, 
they are used in approximation theory as radial basis functions for interpolating scattered 
data on spherical domains, see for example Xu and Cheney, 1992,Fasshauer and Schumaker, 
1998 or Cavoretto and De Rossi, 2010. 

Recently, Gneiting Gneiting, 2013 has reviewed conditions for functions to belong to *&d 
or and uses them to study parametric families of isotropic and stationary correlation 
functions on the sphere. In his work, he also states several problems for future research, 
one of which has been solved in Ziegel, 2013, and the solution to another one is given here. 

The following important characterization of the members of ^ and (due to Schoen- 
berg Schoenberg, 1942 and Chen Chen, Menegatto, and Sun, 2003) plays a key role. It says 
that the members of ^ can be expressed as series expansions using Gegenbauer polynomials 
C^- In particular, the class d > 1, consists of the functions of the form 

♦(•l^^'CiS 9 ^ IM, (2) 

n=0 UJ 

with b n ^ > and Y^=o^n,d = 1- For d > 2 the class Vt^ consists of functions of ^ 
where b U) d > for infinitely many even and infinitely many odd integers n. Consequently, 
it is possible to study the members of ^ or ^d v * a ^ ne coefficients b n ^ of the Gegenbauer 
expansion (2) which are called the d- dimensional Schoenberg coefficients of if). Since 

$lDl 2 Df 3 D..., 

every function in ^d, d > 1, allows a Gegenbauer expansion in terms of C^(cos#) = 
cos(n(9), n = 0, 1, 2, . . . and b Ut i, which is a Fourier cosine expansion. Similarly, if d > 2 

1/2 

every function if) G ^d allows an expansion in terms of Legendre polynomials, since Cn — 
P n is a Legendre polynomial. Hence, it is interesting to ask if it is possible to express higher 
dimensional Schoenberg coefficients in terms of Fourier or Legendre coefficients. 

Referring to Wendland, 1995 we will call the expression of higher dimensional Schoenberg 
coefficients in terms of lower dimensional ones a dimension walk, because this describes a 
method which connects Gegenbauer expansions for different dimensions. In contrast, Wend- 
land uses this phrase to describe the role of the operators Montee and Descente introduced 
by Matheron Matheron, 1965, which work in different ways (see Daley and Porcu, 2013 for 
a characterization of these operators in terms of Schoenberg measures) . 

In general, dimension walks are interesting for studying properties of functions, because 
they allow transferring properties from a certain dimension to another. 

The problem of expressing even and odd dimensional Schoenberg coefficients in terms of 
Fourier and Legendre coefficients can be answered using the following recursive identities, 
stated as Corollary 3 in Gneiting, 2013, which provide a connection between d-dimensional 
Schoenberg coefficients and lower dimensional ones. In particular, for all integers n > 1 it 
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is true that 

bo,3 = b ,i - -b 2 ,i and o„ >3 = -(n + l)(6 n ,i - b n+2 ,i). (3) 

Furthermore, if d > 2, then for all integers n > 

_ (n + d-l)(n + d) (n + l)(n + 2) 

6 ™' d+2 " d(2n + d-l) &n ' d " d(2n + d + 3) K+24 (4) 

These recursive relationships suggests that it is possible to express b n ^k+i,k > 1, as a 
linear combination of Fourier coefficients o n ,i 5 &n+2,i 5 • • • , 0n+2fc,i- Similarly, we can express 
b n ,2k+2, k > 1, as a linear combination of Legendre coefficients b n ^, 0n+2,2> ■ ■ ■ > »n+2fc,2- 

The aim of this work is to provide closed form expressions of the coefficients appearing 
in these linear combinations, which was stated as Problem 2 in Gneiting, 2013. 

2 Main results 

In this section we give the explicit linear combinations for Gegenbauer coefficients in terms 
of Fourier cosine and Legendre coefficients. The proofs can be found in Sections 4 and 5. 

Theorem 2.1. For integers k > 1 andn > the Gegenbauer coefficient b Ht 2k+i of a function 
ip G ^2k+i can be expressed in terms of its Fourier cosine coefficients b n> i, 6 n +2,i> • • • , 0n+2fc.i> 
in that 



b n ,2k+i = y^Qt(" 7 k)b n+ 2j, 



i=0 



where the coefficients a,i(n,k) are given by 

^*)=^(*) (B r* )(n 1 ; M)( r + '^- 1 ' km^m- m 

2 k VV (2A-1)!! (n + *)( fc+ i) (5) 
ao(0) fe) = 1 /or i = n = 0, 

where (2k— 1)!! = n!=i(2* — 1) anc ^ ( x )(m) = aj(sc+l) • • • (x+m — 1) denotes the Pochhammer 
symbol. 

Example 2.2. Consider k = 4. For n > we get 

o (n,4) = «(n + 4)(n + 5)(n + 6)(n + 7) 
ai (n, 4) = -4/c(n + 2)(n + 4)(n + 6)(n + 7) 
o 2 (n, 4) = 6«(n + l)(n + 4) 2 (re + 7) 
a 3 (n,4) = -Ak(u + l)(n + 2)(n + 4)(n + 6), 
04 (n, 4) = re(n+ l)(n + 2)(n + 3)(n + 4), 

where k = ^ 
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One sees that ao(n.,4) and 04(71, 4) can be expressed in simpler forms. In general, for 
i = and i = k equation (5) reduces to 



a*(n,*)=(~) (2i^i)!T (n + 1) W- 

It might be interesting to note that the value of X^i=o °«( n ; i s either or |, despite the 
cumbersome appearance of the explicit formula for ctj(n, k). 

Proposition 2.3. For all integers k > it is 

i=0 I 2 ' " - U 

The proof can be found in Section 4. 

Now let us turn to the analogous problem of finding a linear combination for b n ^k+2i k > 1, 
in terms of Legendre coefficients 6 nj 2, • • ■ , &n+2fc,2- This is very similar to the problem above 
and one suspects similar coefficients. 

Theorem 2.4. For all integers k > 1 and n > we have 

k 

K,2k+2 = E U i( n ' k)b n+2i ,2, (6) 
i=0 

where the coefficients Ui(n,k) are given by 

Ul{n ' ^ = ^^W^O \n ") {n + i + l/2)( fc _i)(n + k + 3/2)^ ' (?) 



3 Applications 

3.1 Continuity 

As mentioned in the introduction, several results are known that connect the behaviour 
of Legendre and Fourier expansion coefficients with the local behaviour of the correlation 
function. Boas Boas, 1967 analyzed Fourier cosine series with non-negative coefficients, 
which correspond to members of If b n 1 4 0, he showed that 

(i) tp € ^1 is Holder continuous with exponent 7 € (0, 1) if and only if b nt i = 0(n _7_1 ) 
and 

(ii) ip has fractal index 7 E (0, 1), meaning that 1 — i/j(0) = 0{9 1 ) as 9 \. 0, if and only if 
b nA = Oin^" 1 ). 

These results allow the construction of members of ^fc+ij k > 0, with a given fractal index 
7 € (0, 1). Note that this construction principle works only if the constructed coefficients 
b n ,2k+i are again non-negative, which might be hard to verify. 
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Example 3.1. Let 

6 

b n ,i = —9-0, n > 

and 60,1 = 0. Obviously, the corresponding function is in ^f. For n,/c > la symbolical 
calculation with Mathematica yields 

" 3k(n + k)B(n/2,k) 2 

b n> 2k + i = X, ai (n,k)b n+2i ,i - n7r 2 {n + 2k) 2 B{ ^ 2k y 

where B(x,y) denotes the Beta function. Note that fo n ,2fc+i > for all n,k > 1 showing 
that the corresponding function is in . 

By Stirling's formula for Gamma functions (see 6.1.37 in Abramowitz and Stegun, 1972), 
we get for fixed y that 

yy/2-y-x ^ 1 + Q( X ~ 1 ) 



B(x,y) =r(y)e» 1 + ^ x 



x) l + Odx + y)- 1 ) 

Since (l + -) 1//2 y x e~ v as x — > 00 we see that 

B(x,y)^T(y)x-y, 

for x — > 00, yielding immediately 

£(n/2,£;) 2 



B(n,2k) 



0{1). 



Consequently, b n ^k+i = 0{n 2 ) for all /c > 1, implying that the higher-dimensional Schoen- 
berg coefficients show the same asymptotic behaviour as b n \. 



3.2 Differentiability 

Hitczenko and Stein Hitczenko and Stein, 2012 derive for a certain class of stationary and 
isotropic random fields on the sphere § 2 necessary and sufficient conditions for their corre- 
sponding two-dimensional Schoenberg coefficients. In particular, they considered the class 
of stationary and isotropic random fields on S 2 with correlation function 

00 9-1-1 

i>(0) = Y,c(n)^-P n (co S 9), (8) 

n=0 

such that the coefficients c(n) fulfill the following: there exist c n > 0, n > 0, such that for 
n > we have c(n) = ^i 2 ^, e > 0, and for n = it is c(0) = Co; furthermore c n is bounded 
below (above) by Ai > (A2 > 0) and the limit linin^oo c n = c exists. They showed that the 
corresponding isotropic and stationary process is A:-times mean square differentiable along 
any fixed direction if and only if e > 2k. 

Their proof is based on showing that the corresponding isotropic covariance function is 
2/c-times differentiable at zero. We can use this result to construct stochastic processes on 
spheres in even dimensions which have k partial derivatives in any direction. Since we know 
that 6 n> 2 = c(n) 2 ^^- we can transfer the assumptions on c(n) readily to b n ^ and with this 
we get a construction principle for higher dimensions. 



5 



4 Proof of Theorem 2.1 and Proposition 2.3 

Proof of Theorem 2.1: First we consider the case n > 0. We proceed by induction for 
k > 1. 

Let /e = 1. For all n > 1 we have 

K,z = 7^(n + l)0n,i - K+2,l), 

yielding ao(n, 1) = ^(n + 1) and a\(n, 1) = — ^(n + 1). Inserting z = and k = 1 into 
formula (5) immediately yields the same results, proving the claim for k = 1. 

Suppose we have proven (5) for an arbitrary k > 1. From this we want to deduce (5) for 
k + 1. With (4) we see by comparing coefficients that 

(n + 2k) {2) + l)( 2 ) , 

b «. 2 (*+i)+i " 2(2k + l)(n + k) bn > 2h+1 " 2(2Hl)(n + H2) 6 " +W 

(n + 2fc) (2) ^ (n + l) (2) ^ 

= 2(2fc + l)(n + fc) ^ ^ )6 " +2i ' 1 " 2(2fc + l)(n + fc + 2) ^ + 2 ' 

i=0 i=0 

(n + 2fc) (2) , ,v, (w + l) (2) , ^ . 

= 2(2fc + l)(n + fc) ao(n ' ^ " 2(2fe + l)(n + /c + 2) afc(n + 2 ' fc ) W+D,i 



+ 2(2fc + 1) ^ 



— ai(n, fc) — — — ai_i(n + 2, fe) 

n + k n + k + I 



i=i 

Using the induction hypothesis and the trivial identity 

+ k)(i) = for integers k, I > 0, (9) 

we see that 

(n + 2fc) (2 ) (n + 2/c) (2 ) 1 (n + fc)n(n + l)( 2fe _i) 

ao{n, K) 



2(2k + l)(n + k) uv ' ; 2(2k + l)(n + k) 2 k (2k - l)!!(n) (fc+1) 

1 ra(ra + /c + l)(n + l)( 2A . +1 ) 
"2^1 (2fc + l)!!(ra) (fe+2) ' 

proving the validity of (5) for i = and 

- , , (n + 1)(2) , -a fc (n + 2,A;) 

2(2k + l)(n + k + 2) ky ' ; 

(n + 1) (2 ) (-l) fc (n + 2 + fc)(n + 2 + 2fc)(n + 2 + l) (2fc _i) 

2(2fc + l)(n + k + 2) 2 k (2k - l)!!(n + 2 + fc)(fc+i) 

(_l)*+i ( n + 2 + 2fc)(n + l) (2fc+1) 



2*+* (2fc + l)!!(n + 2 + % +1) 
(_l)fc+i („ + fc + l)(n + 2k + 2)(n + l) (2fc+1) 
2W (2fc + l) \l(n + l + % +2) 

which shows the validity of (5) for z = k + 1. 



It remains to show for 1 < % < k that 

(-1)* (k + 1\ (n + k + l)(n + 2i)(n + l) (2 fc+i) 



2^+1 V i J (2fc + l)!!(n + i) (fc+2) 



1 



2(2/c + 1) 



( n + 2fc )(2) , (" + %) , . 

— ai(n, k) v > i_i(n + 2, fc) 

n+K n+K+2 



where dj(n, fc) and Oj_i(n + 2, A;) can be expressed as in (5). Plugging the induction hypoth- 
esis into (10) we can reformulate this as 

{-I) 1 (k + 1\ (n + k + l)(n + 2i)(n + l) (2fc+ i) 



V i y (2A; + l)!!(n + z) (fc+2) 
_ (-!)< (n + 2k) {2) (k\ (n + k)(n + 2i){n + l) {2 k-i) 



2 k + 1 n + k \ij (2fc + l)!!(n + i) (fe+1) 
(_l)*-i (n + l) (2) / k \(n + 2 + k)(n + 2i)(n + 3) {2k _ 1) 



2 k + 1 n + k + 2\i-lj (2k + l)!!(n + 1 + *)(fc+i) 
By using (9) and canceling factors we see that this is equivalent to 

k + l\. ^ 1 AA 1 f k \ 

(n + k + 1)- = . + 



i 



After multiplying with (n + i)(fc+ 2 ) it remains to show that 

k+ i i yn + k + l)=Q(n + l + k + l)+( i k _ i 

The right hand side of (11) equals 

k\ /n + i + k+1 n + i 

+ 



(k — — 1)! V i k-i + 1 

k\ 



(k + l-i)\i\ 



[(k - % + l)(n + % + k + 1) + i(n + i) 



(k + l)(n + k + 1) = f * + (n + k + 1), 



showing the validity of (11) and we are done. 

Now consider the case n = 0. As (3) and (4) show, the coefficient 6 n d+2 satisfies for all 
n > and all d > 2 the same recursive property. This shows that aj(0, fc) is given by (5) for 
all z > 0. If n = we can reformulate (4) and together with (3) we see that the following 
holds for all d > 1 



?0,d+2 = 00,d 



d(d + 3) z ' a ' 

Using this relationship recursively, we see that for all d > 1 it holds that 

bo,d+2 = &o,i — R} 

where the remainder term R does not depend on &o,i- This shows that ao(0,A;) = 1 for all 
k > 1. □ 
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Proof of Proposition 2.3: Let n > 0. Using (n + i)(k+i) = ( ) (k + l)\ and Theorem 
2.1 we see that 



T «*(», fc) = (n + "l ( : + 1 ! ( ,r i) E (-ir f fc ) ? n \ 2t 

h 2 fe (2fc-l)!! ' W(" + *)(*- 



i=0 v ' i=0 



V (« + %+!) 



(n + fc)(n + l)( 2 fc-i) v^, ifk\ n + 2i 



' 2*(fc + l)!(2fc-l)!! ^ ; W 
Hence, it suffices to prove that 

yVi)^ ra + 2i =o 

i=o v 7 \ fc+l / 

which is equivalent to 

k 



I fc+l / i=0 V 7 V fc+l / 



i=0 

Now the left-hand side of (12) equals 

^Wr+tf) 2fc+i ( ^ r) ( ^)' 

where the first equality is due to the following result of R. Frisch which can be found, for 
example, as Note 21 in Netto, 1927 



i=0 x%/ \ c 



h+ c) k + C (b-c) ' 



where b > c are positive integers. 

For the right-hand side of (12) we get in a very similar way 



/n+i+fcN / (Jb-i)!(*-l)! f, 

i=0 v 7 V fc+l J i=0 y ' V ' V fc+l 

= -2* £(-!)< 



'fc-l\ 1 



i=0 
fc 



7 I fc+l J 
1 



fn+«+l+fc\ 
I fc+l / 



nl k + \ 1 fc + l 

2fc- 



2/j (2k+n^ ^2fc+n^ 
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where we use (13) for the final equality, thereby showing (12). 
Now consider the case n = 0. For i > 0, equation (5) simplifies to 



<n{0,k) = (-1)' 



k\ 1 



which is also valid for i = 0, since in this case it reduces to 1. Hence, by using (13) with 
b = c = k we get 

k k 



S „ i( o,*) = B-i) i Q) 7E k = i- 

i=0 i=0 v 7 V k I 



□ 



5 Proof of Theorem 2.4 

Proof of Theorem 2.4: We proceed by induction. Let k = 1. Then 



h n,i = -An + 1)( 2 ) 



1 . 1 



n,2 — n — rOn+2,2 



2n+l ' 2n + 5 



implying tt (n, 1) = 2 (2n+i) ( re + 1 )(2) and x ) = ~ 2(2n+5) ( n + 1 )(2)- Inserting A; = 1, 

and i = and i = 1 respectively in equation (7) proves the claim for k = 1. 

Suppose we have proven (7) for a fc > 1. We use this to show the validity of (7) for k + 1. 
Using (3) and (4) and the induction hypothesis we see that 



n 



+ 2fc + l) (2) (« + l)(2) 



&n,2( fc+ i)+2 - 2 (fe + l)(2n + 2fc + l) &n ' 2fc+2 " 2(fc + l)(2n + 2A ; + 5) 0n+2 ' 2fc+2 

(n + 2fc + l) (2) ^ (n+l) (2) * 

= 2(fc + l)(2n + 2A: + l) ^ ^ ^ +2 *' 2 " 2(* + l)(2n + 2fc + 5) ^ ^ + fc ^+ 2 + 2 *' 2 

(n + 2fc + l) (2 ) ( n + !)(2) 

-u (n, fc)6 n , 2 - , lV ~ , „, , ^ "fcfr* + 2, fc)&n + 2(fc+i),2 



2(A: + l)(2n + 2fe + 1) v ' ; ' 2(fc + l)(2n + 2A; + 5) 



i=l 



^ + ^ + 1)(2) "«<(», k) - „„ , , K ^ _ l( n + 2, fc) 



2(fc + l)(2n + 2A; + 1) v ' ' 2(k + l)(2n + 2A; + 5) 



&n+2i, 



2- 



Using (2/c — 1)!! = and (9) we see that 



(TO + 2fc + l) (2) _ (ra + 2fc + l) (2) (2fc-l)!! /2fc + n\ 1 

2(k + l)(2n + 2k + l) U ^ n ' ' ~ 2 k + 1 (k + l)(2n + 2k + 1) V n 7(n + l/2) (fc) 



1 (2k + 2 + n)l 1 
2fe+i2fc+i ( k + l)\n\ (n + l/2) (fc+1) 

(2fc + l)!! / / 2£; + 2 + n\ 1 



2 fc+1 V n y (n + l/2) (fc+1) ' 
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proving the claim for % = 0. In a very similar way we get 
(n + 1)(2) 
~2{k + l){2n + 2k + h) Uk{n + 2 ' k) 
_ , )k+1 (ra + l) (2) (2fc-l)!! /2k + n + 2\ 
( ' 2 fe + 1 (ft + l)(2?i + 2ft + 5) \ n + 2 ) 
= ( 1 {2k + 2 + n)\ 1 



in + ft + 7/2) {k) 



2 k+i 2 k+i (k + l)\n\ (n + fc + 5/2)( fc / 2 ) 



fc+1 (2ft + 1)!! /2ft + 2 + n^ 1 



2 fc+1 V " /(n + fc + 5/2) (fc/2) ' 



confirming the claim for i = k+l. 

Now let 1 < i < k. We have to show that 

. _(2fc + l)!! ^ + lV2H2 + n 

iti(n,fc + 1) = (-1) 4 ' 



2 fc+1 V » A « / (n + t + l/2) (fc+1 _ (n + fc + l + 3/2) (i) 

( n + 2 fc + l) {2) (n + 1) ( 

Ui(n, k) - — — — — — Ui-i(n + 2, ft), 



2(k + l)(2n + 2k + 1) v ' ' 2(k + l)(2n + 2fc + 5) 
which is equivalent to 
;(2ft + l)!! fk + l\ 



2 fc+1 V i J\ n J (n + i + l/2) (fc+1 _ i) (n + fc + 1 + 3/2) (i) 
(-l)Hn + 2ft + l) (2 )(2fc - 1)!! /ft\ /2ft + n\ 1 

2 fc + 1 (ft + l)(2n + 2ft + l) V*/ V n ) (n + i + l/2) (ft _ (n + A; + 3/2) (i) 
(-l) i - 1 (n + l) (2) (2fc-l)!! / ft \/2fc + n + 2\ 1 



2 fc + x (fc + l)(2ra + 2fc + 5) V*-VV n + 2 / ( n + * + 1 + V%-i+i) ( n + 2 + ft + 3/2) (i _ 1} 

(14) 

Using (9) and doing similar manipulations as in the cases i = and i = ft + 1 we see that 
(14) is equivalent to 

k + l\ 1 



i ){n + i + 1/2)^+!^) (n + k + 5/2) (i) 
ft\ 1 / 



^ (n + t + l/2) (fc+1 _ (n + ft + 3/2) (i) \i - 1) (n + i + 3/2) (fc+1 _ i} (n + ft + 5/2) (i) 
Multiplying with (n + i + l/2)^ +1 _j)(n + ft + 5/ 2 )(i) illustrates that we have to show 
k + l\ _ fk\n+k + i + 3/2 ^ ( k \n + i + 1/2 



i J \ij n + k + 3/2 \i - lj n + k + 3/2 ' 
Simplifying the right hand side yields 

'k\(n + k + i + 3/2) (Jfc - i + 1) + i(n + i + 1/2) 
(n + ft + 3/2)(fc-i + l) 

'Jfe\ (ft+ l)(n + ft + 3/2) /ft + 1 



(n + ft + 3/2)(ft-i + l) V * 
and we are done. □ 
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